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Abstract. 

With the advent of galaxy surveys which provide large samples of galaxies or 
galaxy clusters over a volume comparable to the horizon size (SDSS-III, 
HETDEX, Euclid, IDEM, LSST, Pan-STARRS, CIP etc.) or mass-selected large 
cluster samples over a large fraction of the extra-galactic sky (Planck, SPT, ACT, 
CMBPol, B-Pol), it is timely to investigate what constraints these surveys can 
impose on primordial non-Gaussianity. I illustrate here three different 
approaches: higher-order correlations of the three dimensional galaxy 
distribution, abundance of rare objects (extrema of the density distribution), and 
the large-scale clustering of halos (peaks of the density distribution). Each of 
these avenues has its own advantages, but, more importantly, these approaches 
are highly complementary under many respects. 
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1 Introduction 

The recent advances in the understanding of 
the origin and evolution of the Universe have 
been driven by the advent of high-quality data, 
in unprecedented amount (just think of WMAP 
and SDSS for example). Despite this, most 
of the information about cosmological param- 
eters come from the analysis of a massive com- 
pression of the data: the power spectrum of 
their statistical fluctuations over the mean. The 
power spectrum is a complete statistical de- 
scription of a random field only if it is Gaus- 
sian. 

Even the simplest inflationary models predict 
deviations from Gaussian initial conditions. 
These deviations are expected to be small, al- 
though "small" in some models may be "de- 
tectable". For a thorough review of inflationary 
non-Gaussianity see [IJ, for our purpose it will 
be sufficient to say that to describe inflation- 
motivated departures from gaussian initial con- 
ditions many write (ESUllIl: 

O = + /;Vl(0'-(0')) . (1) 

Here (p denotes a gaussian field and O 
denotes Bardeen's gauge-invariant potential, 
which, on sub-Hubble scales reduces to the 
usual Newtonian peculiar gravitational poten- 
tial, up to a minus sign. In the literature, there 
are two conventions for Eq. ([I]): the large-scale 
structure (LSS) and the Cosmic Microwave 
Background (CMB) one. In the LSS conven- 
tion O is linearly extrapolated at z = 0; in the 
CMB convention O is instead primordial: thus 
fil' = = oo)/g(0)/j^f^ ~ 1.3/^f^ where 
g{z) denotes the linear growth suppression fac- 
tor relative to an Einstein-de-Sitter Universes. 
In the past few years it has become customary 
to always report f^^^ values even if, for sim- 
plicity as it will be clear below, one carries out 
the calculations with f^l^. 

While for simplicity one may just assume /ml 
in eq. [T] to be a constant (yielding the so-called 
local model or local-type) in reality the expres- 
sion is more complicated and /nl is scale and 



configuration dependent. In general, the non- 
Gaussianity is specified by writing down the 
bispectrum of O. For example, one can see that 
for the local model the bispectrum is: 

B^(kuk2,ks) = 2fni^P^(ki)P^(k2) + 2cyc. (2) 

Where P denotes the power spectrum and it is 
often assumed that P^ = Pd,; "eye." denotes two 
cyclic terms over ki,k2, k^. 

It has been shown [6] that for non- 
Gaussianity of the local type the bispectrum is 
dominated by the so-called squeezed configu- 
rations, triangles where one wavevector length 
is much smaller than the other two. Models 
such as the curvaton for example have a non- 
Gaussianity of the local type. Standard, single- 
field slow roll inflation also yields a local non- 
Gaussianity but with an unmeasurably small 
/nl (see [IJ and references therein). On the 
other hand, many inflationary models have an 
equilateral-type non Gaussianity, i.e. the bis- 
pectrum is dominated by equilateral triangles 
[6]. Ref. [[7l[8l have proposed a functional form 
which closely approximates the behavior of the 
inflationary bispectrum and which is useful for 
efficient data- analysis: 

BikuhX) = 6fnL{-P(h)P(k2) + 2cyc. 

- 2[P(h)P(k2)P(ks)f' (3) 
+ P''\h)P^I\k2)P{h) + 5cyc.) . 

Note that the same numerical value for /nl 
gives rise to a larger skewness in the local case 
than in the equilateral case (e.g. see [9|), ex- 
plaining why the CMB constraints on /nl are 
weaker for the equilateral case iTTllSlfTOll. 

Specific deviations from a single field, slow 
roll, canonical kinetic energy, Bunch-Davies 
vacuum, leave their specific signature on the 
bispectrum "shape" (i.e. the dependence of B 
on the shape of the triangle made by the three k 
vectors), see discussion in IfTTI and references 
therein. 

Non-Gaussianity therefore offers a probe of 
aspects of inflation (namely the interactions 



of the inflaton) that are difficuk to probe by 
other means (i.e., measuring the shape of the 
primordial power spectrum and properties of 
the stochastic background of gravity waves). 
So, how could primordial non-Gaussianity be 
tested? 

One could look at the early Universe: by 
looking at CMB anisotropics we can probe cos- 
mic fluctuations at a time when their statisti- 
cal distribution should have been close to their 
original form but the signal is small. On the 
other hand, one could analyze the statistics of 
the large-scale structures, close to the present- 
day, when the signal is larger, but this is a more 
complicated approach, since gravitational in- 
stability (for the dark matter distribution) and 
bias (for galaxies or clusters of galaxies) in- 
troduce non-Gaussian features in an initially 
Gaussian field and they mask the signal one 
is after. Finally, the abundance of rare events 
(such as galaxy clusters and high-reshift galax- 
ies) probes the tails of the PDF of the density 
field, which are extremely sensitives to devia- 
tions from Gaussianity. Here, I will concentrate 
on the signature of non-Gaussianity on large- 
scale structure (i.e. at redshift z < 1) as they 
can be traced by galaxy surveys (i.e., I will not 
consider wide field weak gravitational lensing 
surveys); other contributions to this review will 
focus on non-Gaussian signatures on the CMB, 
thus here it will be sufficient to give only a very 
brief and succinct introduction to the subject. 

At recombination the density fluctuations are 
small( 6^ ~ 10"^), and the CMB temperature 
fluctuations are directly related to O making 
this a very clean probe. However, effectively 
only one redshift can be tested giving us only a 
2-dimensional information. 

The most widespread technique for testing 
Gaussianity in the CMB is to use the CMB bis- 
pectrum: 

/ m, m2 _ p (Ci f 3 \ (A\ 

mi mi ni-i J \ 

where the a'" are the coefficients of the spher- 
ical harmonic expansion of the CMB temper- 



ature fluctuation: AT/T = 2^^^/"^"' ^iid the 
presence of the 3- J symbol ensures that the 
bispectrum is defined if li + I2 + h = even, 
+ 4 > li > \lj - €k\ (triangle rule) and 
that mi + m2 + = 0. It should be how- 
ever clear that secondary CMB anisotropics 
and foregrounds also induce a CMB bispectrum 
which can mask or partially mimic the signal 
see e.g., Refs. [HI El [H [H [Bl [H [HI and 
references therein. 

In the last few years, this area of research 
has received an impulse, motivated by the re- 
cent full sky CMB data from WMAP. In par- 
ticular it has been shown that the constraints 
can be greatly improved by effectively "recon- 
structing" the potential O from CMB tempera- 
ture and polarization data rather than simply us- 
ing the temperature bispectrum alone ilEl IT91 . 
This technique would yield constraints on non- 
Gaussianity of the local type of A/nl ~ 1 for an 
ideal experiment and A/nl ~ 3 for the Planck 
satellite. This is particularly promising as /nl 
of order unity or larger is produced by broad 
classes of inflationary models (see e.g., p"J and 
references therein). 

Currently, the most stringent constraints for 
the local type are 27 < /nl < 147 at the 95% 
confidence (central value 87) from WMAP 3 yr 
data [EOl; and from the WMAP 5 years data, 
-9 < /nl < 111 at the 95% confidence level 
(central value 55) [ 10| and -4 < /nl < 80 [21J. 
Despite the heated debate on wether /nl = is 
ruled out or not, the two measurement are not 
necessarily in conflict: the two central values 
differ by only about Icr; different, although not 
independent, data sets were used with different 
galactic cuts, and the maximum multipole con- 
sidered in the analyses is also different. What 
makes the subject very interesting, is that, if the 
central value for fyi is truly around 60, forth- 
coming data will yield a highly- significant de- 
tection. 



2 Higher-order correlations 

Theoretical considerations (see discussion in 
flllil and references therein) lead us to de- 
fine primordial non-Gaussianity by its bispec- 
trum. While in principle there may be types of 
non-Gaussianity which would be more directly 
related to higher-order correlations (e.g., [[22l| 
and references therein), and while a full de- 
scription of a non-Gaussian distribution would 
require the specification of all the higher-order 
correlations, it is clear that quantities such as 
the bispectrum enclose information about the 
phase correlation between fc-modes. In the 
Gaussian case, different Fourier mode are un- 
corrected (by definition of Gaussian random 
phases) and a statistic like the power spectrum 
does not carry information about phases. The 
bispectrum is the lowest-order correlation with 
zero expectation value in a Gaussian random 
field. But, even if the initial conditions were 
Gaussian, non-linear evolution due to gravita- 
tional instability generates a non-zero bispec- 
trum. In particular, gravitational instability has 
its own "signature" bispectrum, at least in the 
next-to-leading order in cosmological perturba- 
tion theory nil: 

B{k,XX) = 2P{h),P{k2)J{hX) + 2cyc. 

(5) 

where J{ki,k2) is the gravitational instability 
"kernel" which depends very weakly on cos- 
mology and for an Einstein-de-Sitter Universe 
is: 



^ ^ 5 k\ -ki 
Jikuk2) = - + ^ 
7 2kik2 
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(6) 

In the highly non-linear regime the detailed 
form of the kernel changes, but it is something 
that could be computed and calibrated by ex- 
tending perturbation theory beyond the next-to- 
leading order and by comparing with numer- 
ical N-body simulations (see other contribu- 
tions in this issue). It was recognized a decade 
ago L4J that this signal is quite large compared 



to any expected primordial non-Gaussianity 
and that the primordial signal "redshifts away" 
compared to the gravitational signal. In fact, 
a primordial signal given by a local type of 
non-Gaussianity parameterized by a given /nl, 
would afi'ect the late-time dark matter density 
bispectrum with a contribution of the form 



Bf^^i'"'-"'{kuk2%,7) 



(7) 



2fnLP(h)P(k2)- 
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+ 2cyc. 



Z)(z)/D(z = 0) 
where D(z) is the linear growth function which 
in an Einstein-de Sitter universe goes like (1 -l- 
z)~^ and 



Mih) 
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T(k) denoting the transfer function, Hq the 
Hubble parameter and Q.^^) the matter density 
parameter. Clearly the two contributions have 
different scale and redshift dependence and the 
two kernel shapes in configuration space are 
different, thus, making the two component, at 
least in principle and for high signal-to-noise, 
separable. 

Unfortunately, with galaxy surveys, one does 
not observe the dark matter distribution di- 
rectly. Dark matter halos are believed to be 
hosts for galaxy formation, and different galax- 
ies at different redshifts populate halos fol- 
lowing different prescriptions. In large-scale 
structure studies, often the assumption of lin- 
ear, scale independent bias is made. A lin- 
ear bias will not introduce a non-zero bispec- 
trum in a Gaussian field and its effect on a 
field with a non-zero bispectrum is only to re- 
scale its bispectrum amplitude. This is, how- 
ever, an approximation, possibly roughly valid 
at large scales for dark matter halos, and when 
looking at the power spectrum, but unlikely 
to be true in detail. To go beyond the lin- 
ear bias assumption, often the assumption of 
quadratic bias is made, where the relation be- 
tween dark matter overdensity field and galaxy 



field is specified by two parameters: bi and 

and b2 are assumed to be scale-independent 
(although this assumption must break down 
at some point) but they can vary with red- 
shift. Clearly, a quadratic bias will introduce 
non-Gaussianity even on an initially Gaussian 
field. In summary, for local non-Gaussianity 
and scale-independent quadratic bias we have 

mm-. 

B(lluL L z) = 2Piki)Pik2)biiz)' X (9) 

ni^iJi) , b2(z)] 

-^NL p., , + J(ki,k2) + T— — +cyc. 
D(z) 2oi(z) 

Before the above expression can be compared 
to observations it needs to be further compli- 
cated by redshift space distortions (and shot 
noise). Realistic surveys use the redshift as a 
proxy for distance, but gravitationally-induced 
peculiar velocities distort the redshift-space 
galaxy distribution. We will not go into these 
details here as including redshift space distor- 
tions (and shot noise) will not change the gist 
of the message. 

From a practical point of view, it is important 
to note that photometric surveys, although in 
general can cover larger volumes that spectro- 
scopic ones, are not suited for this analysis: the 
projection effects due to the photo-z smearing 
along the line-of-sight is expected to suppress 
significantly the sensitivity of the measured bis- 
pectrum to the shape of the primordial one (see 
e.g.,[|251|26l). 

Ref. H concluded that "CMB is likely to 
provide a better probe of such [local] non- 
Gaussianity". Much more recently, ETl revis- 
ited the issue and found that, assuming a given 
-known- redshift dependence of the {bi, b2) bias 
parameters and an all sky survey from z = to 
z = 5 with a galaxy number density of at least 
5 X 10~'*/z^/Mpc^, the galaxy bispectrum can 
provide constraints on the /nl parameter com- 
petitive with CMB. However, for all planned 
surveys, the forecasted errors are much larger 



than Planck forecasted errors. This holds qual- 
itatively also for the equilateral case. 

While the gravitationally-induced non- 
Gaussian signal in the bispectrum has been 
detected to high statistical significance (see 
[|28l and references therein; other contributions 
to this issue); the non-linear bias signature 
is not uncontroversial, and there have been 
so far no detection of any extra (primordial) 
bispectrum contributions. 

Of course one could also consider higher- 
order correlations. One of the advantages of 
considering e.g., the trispectrum is that, con- 
trary to the bispectrum, it has very weak non- 
linear growth [29], but has the disadvantage that 
the signal is de-localized: the number of possi- 
ble configurations grows fast with the dimen- 
sionality n of the n-point function! 

In summary, higher-order correlations as ob- 
served in the CMB or in the evolved Universe, 
can be used to determine the bispectrum shape. 
The two approaches should be seen as comple- 
mentary as they are aff"ected by diff"erent sys- 
tematic eff"ects and probe diff"erent scales. The 
next two probes we consider have a less rich 
sensitivity to the bispectrum shape, but their 
own peculiar advantages. 

3 The Mass Function 

The abundance of collapsed objects (dark mat- 
ter halos as traced e.g., by galaxies and galaxy 
clusters) contains important information about 
the properties of initial conditions on galaxy 
and clusters scales. The Gaussian assumption 
plays a central role in analytical predictions for 
the abundance and statistical properties of the 
first objects to collapse in the Universe. In this 
context, the formalism proposed by Press & 
Schechter [30], with its later extensions and im- 
provements, has become the standard lore for 
predicting the number of collapsed dark matter 
halos as a function of redshift. However, even 
a small deviation from Gaussianity can have a 
deep impact on those statistics which probe the 
tails of the distribution. This is indeed the case 



for the abundance of high-redshift objects like 
galaxies and clusters at z > 1 which correspond 
to high peaks, i.e. rare events, in the underly- 
ing dark matter density field. Therefore, even 
small deviations from Gaussianity might be po- 
tentially detectable by looking at the statistics 
of high-redshift systems. Before proceeding let 
us introduce some definitions. 

We are interested in predictions for rare ob- 
jects, that is the collapsed objects that form in 
extreme peaks of the density field 6(x) = 6p/p. 
The statistics of collapsed objects can be de- 
scribed by the statistics of the density pertur- 
bation smoothed on some length scale R (or 
equivalently a mass scale M = 4/3nR^p), 6r. 

To incorporate non-Gaussian initial condi- 
tions into predictions for the smoothed density 
field, we need an expression for the probability 
distribution function (PDF) for 6r. For a partic- 
ular real-space expansion like Eq. ([T]), one may 
make a formal change of variable in the Gaus- 
sian PDF to generate a normalized distribution 
on . However, this may not be possible in gen- 
eral and the change of variables does not work 
for the smoothed cumulants of the density field. 
In general, the PDF for a generic non-Gaussian 
distribution can be written exactly as a function 
of the cumulant's generating function "IVr for 
the smoothed density field: 

(10) 



with 



^RW = y^-^Mn,R (11) 
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Fig. 1 Skewness St,^ of the density field at 
z = as a function of the smoothing scale R 
for different types of non-Gaussianity. Figure 
reproduced from ||32l . 

jx^^R is related to the underlying bispectrum by: 



d\i(Pk2(Pk-i -> -> o 
— r- Bs^R{ki , k2, h)6^ , f 



(12) 

where 6^ denotes the Dirac delta function and 
B^R denotes the bispectrum of the 5 overden- 
sity field smoothed on scale R. It is related 
to the potential one trivially by remembering 
the Poisson equation: 5R(k) = M{k)WR{k)(^(k). 
Here, WR{k) denotes the smoothing kernel: usu- 
ally taken to be the Fourier transform of the top 
hat window. In any practical application the di- 
mensionality of the integration can be reduced 
by collapsing the expressing kj as a function of 
ki and ^2- 

It is important at this point to make a small 
digression to specify definitions of key quan- 
tities. Even in linear theory, the normalized 
skewness of the density field depends on red- 
shift; however in the Press & Schechter frame- 
work one should always use linearly extrapo- 
lated quantities at z = 0. In this context there- 



where ^„r denote the cumulants and for ex- fore, when writing 5 3,^ = /nl'S'^'^^ \if5:('^^ 'is 



'/nl=1 



ample fi2.R = <^r = (^r) ^nd the skewness 
yU3s is related to the normalized skewness of 
the smoothed density field Ss^r = h^^rIix\j^. It 
is useful to define a "skewness per /nl unit" 



'/nl=1 
'3,i? 



SO that 53,R = /nl^" 



/nl=1 
3,R 



The skewness 



that if the density field extrapolated linearly at 
z = then /nl must be the LSS one and not the 
CMB one. 

To compute the abundance of collapsed ob- 
jects fro the PDF one will then follow the Press 



& Schechter swindle: first compute 



n>5Az,R) 



dSi^rm (13) 



where 6c denotes the critical threshold for col- 
lapse; then the number of collapsed objects is: 



n(M, z)dM = 2 



SnGM 



dr{> 6c\z,R) 



dM 



■ (14) 



Note that the redshift dependence is usually en- 
closed only in 6^- (Tm is computed on the field 
linearly extrapolated at z = 0, and (5c(z) = 
Ac(z)D(z = 0)/D(z) and Ac(z) depends very 
weakly on redshift and (z = 0) - 1.68. 

Eq. (fT4l) however cannot be computed an- 
alytically and exactly starting from Eq. (flOl) : 
some approximations need to be done in order 
to obtain an analytically manageable expres- 
sion. Two approaches have been taken so far in 
the literature which we will briefly review be- 
low. 

3.1 MVJ approach. The authors of Ref.[l3ll 
proceed by first performing the integration 
over 6r to obtain an exact expression for 'P{> 
6c\z, M). At this point they expand the generat- 
ing functional to the desired order, e.g., keep- 
ing only terms up to the skewness, then per- 
form a Wick rotation to change variables and 
finally a saddle-point approximation to evalu- 
ate the remaining integral. The saddle point ap- 
proximation is very good for large thresholds 
5cI(Tm » 1, thus for rare and massive peaks. 
For the final expression for the mass function 
they obtain: 



h(M, z) = 2 
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(15) 



where ctm denotes the rms value of the den- 
sity field, the subscript M denotes that the den- 
sity field has been smoothed on a scale R{M) 



corresponding to R{M) = [A/3/(4p)]'^^, and 
= 6c Yi'^-^^5^Im73. 

This derivation shows that the mass function 
in principle depends on all cumulants, but that 
if non-Gaussianity is described by a bispectrum 
(and all higher order connected correlations are 
assumed to be zero or at least negligible), it de- 
pends only on the skewness. The mass func- 
tion does not carry explict information about 
the shape of non-Gaussianity. Nevetheless for a 
given numerical value of /nl the skewness can 
have diff"erent amplitude and scale dependence 
for different models of non-Gaussianity, as il- 
lustrated in Fig. ([U). 

3.2 LMSV approach The authors of Ref. [9] 
(LMSV) instead proceed by using the saddle 
point approximation in the expression for P(,6r) 
and then using the Edgeworth expansion trun- 
cated at the desired order. The resulting sim- 
plified expression for the PDF can then be in- 
tegrated to obtain 'P(> 6c\z,R) and derived to 
obtain the mass function. The final mass func- 
tion in this approximation is given by: 



h(M, z) = 2 
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(16) 



Without knowledge of all higher cumulants one 
is forced to use approximate expressions for the 
mass function. By truncating the Edgeworth 
expansion at the linear order in the skewness, 
the resulting PDF is no longer positive-definite. 
The number of terms that should be kept in the 
expansion depends on 5,.(z)/cr^. The trunca- 
tion used is a good approximation of the true 
PDF if 6c{z)/crM is small (and non-Gaussianity 
is small) but for rare events (the tails of the 
distribution) 6c(z)/crM is large. One thus ex- 
pect this approximation to break down for large 
masses, high redshift and high /nl. 




Fig. 2 Correction to the Gaussian mass function as measured in different non-Gaussian simula- 
tions. There is now agreement between different simulations. The y axis should be interpreted as 
Logio'R{M, z, /nl). Reproduced from fig. 4 and 5 of Ref. [,33J . 



Ref. [9| quantified the range of validity 
of their approximation by assuming that when 
terms proportional to 5 3 become important is 
no longer valid to neglect terms proportional 
to higher-order cumulants. Then they define 
the validity regime of their mass function to 
be where corrections from the S j are unimpor- 
tant. They find, as expected, that for very mas- 
sive objects the approximation breaks down and 
that the upper mass limit for applicability of the 
mass function decreases with redshift and /nl- 
But for low masses, redshifts and /nl their for- 
mula is better than the MVJ. On the other hand 
the MVJ range of validity extends to higher 
masses, redshifts and /nl values, as expected, 
as MVJ applied the saddle point approximation 
to P(> 6c\M, z) which is an increasingly good 
approximation for rare objects. 

Of course, the natural observable to apply 
this method to are not only galaxy surveys (and 
the clusters found there), but, especially suited, 
are the mass-selected large clusters surveys of- 
fered by on-going Sunyaev-Zeldovich experi- 
ments (e.g., Planck, ACT, SPT). 

A detailed comparison with N-body simula- 
tions is the next logical step to pursue. 



3.3 Comparison with N-body simulations Be- 
fore we proceed we should consider that the 
Press & Shechter formulation of the mass func- 
tion even in the Gaussian initial conditions 
case, can be significantly improved see e.g., 
[|34l |35l |36l. Much improved expressions have 
been extensively calibrated on Gaussian initial 
conditions N-body simulations. The major lim- 
itations in both the MVJ and LMSV derivations 
(since they follow the classic Press & Shechter 
formulation) are the assumption of spherical 
collapse and the sharp A;-space filtering. In ad- 
dition, the excursion set improvement on the 
original Press & Schechter swindle relies on the 
random-phase hypothesis, which is not satisfied 
for non-Gaussian initial conditions. Since these 
improvements of the mass function have not yet 
been generalized to generic non-Gaussian ini- 
tial conditions (but work is on-going, see other 
contributions in this issue) the analytical results 
above should be used to model fractional cor- 
rections to the Gaussian case. 

Thus the non-Gaussian mass function, 
nNGiM,z) can be written as a function of a 
Gaussian one, nQ{M,z) (accurately calibrated 
on N-body simulations) with a non-Gaussian 



correction factor H (see e.g., [[37ll9l): 



nNG(M, z) = nciM, z^S^, M, z) (17) 



where 



n(M,z,/NL) 



(18) 



n(M,z,/NL = 0) 

and h is given by the MVJ or LMSV approxi- 
mation. The correction 'R can then be calibrated 
on N-body simulations. 

Ref. [[33ll argue that the same correction 
that in the Gaussian case modifies the collapse 
threshold, 6c, to improve over the original Press 
& Shechter formulation, may apply to the non- 
Gaussian correction. The detailed physical in- 
terpretation of this is still matter of debate in the 
literature Il38l [391 l40l. In summary, Ref. [l33j| 
proposes to write the non-Gaussian correction 
factor for the MVJ (HI case as: 



'Rng(M, z, /nl) = exp 
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where 6cc denotes the modified critical density 
for collapse, which for high peaks is 6cc ~ 
6c -^/q. Ref. [|33l calibrated these expressions 
on N-body simulations to find q = 0.75. We 
anticipate here that the validity of this extrapo- 
lation (i.e. in terms of a correction to the critical 
collapse threshold) can be tested independently 
on the large-scale non-Gaussian halo bias as 
described in ^ Note that, in both cases, in 
the limit of small non-Gaussianty the correction 
factors reduce to 
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Fig. 3 The points show the non-Gaussian cor- 
rection to the mass function as measured in the 
N-body simulations of |[33l . Blue corresponds 
to /nl = 200 and red to /nl = -200. the dashed 
lines correspond to the MVJ formulation and 
the dot-dashed lined to the LMSV formulation. 
In both cases the substitution 6c — > 6ec has 
been performed. The y axis should be inter- 
preted as "RiM, z, /nl = 200). Reproduced from 
fig. 7 Ref. (33. 

Non-Gaussian mass functions have been com- 
puted from simulations and compared with dif- 
ferent theoretical predictions in several works 
[El m |43l m im . In the past, conflicting re- 
sults were reported, but the issue seems to have 
been settled: there is agreement among mass 
function measured from different non-Gaussian 
simulations performed by three different groups 
as shown for example in Fig. (O. As expected 
both MVJ and LMSV prescriptions for the non- 
Gaussian correction to the mass function agree 
with the simulation results, provided one make 
the substitution 6c — > 6ec, with some tenta- 
tive indication that MVJ may be better for very 
massive objects while LMSV performs better 
for less rare events. This is shown in Fig. Q 
where the points represent measurements from 
N-body simulations presented in [[33ll . 

3.4 Voids While galaxy clusters form at the 
highest overdensities of the primordial density 



field and probe the high density tail of the PDF, 
voids form in the low density regions and thus 
probe the low-density tail of the PDF. Most of 
the volume of the evolved universe is under- 
dense, so it seems interesting to pay attention to 
the distribution of underdense regions. A void 
distribution function can be derived in an anal- 
ogous way to the Press Schechter mass func- 
tion by realizing that negative density fluctua- 
tions grow into voids [ [32| : a critical underden- 
sity 6y is necessary for producing a void and 
this plays the role of the critical overdensity 
Sc for producing bound objects (halos). The 
more underdense a void is the more negative 
Sy becomes. The precise value of 6y(z) de- 
pends on the precise definition of a void (and 
may depend on the the observables used to find 
voids); realistic values of 6y(z = 0) are ex- 
pected to be > -1. In the absence of a better 
prescription, here, following [32J, 6y is treated 
as a phenomenological parameter and results 
are shown for a range of 6y values. To derive 
the non-Gaussian void probability function one 
proceeds as above with the only subtlety that 6^ 
is negative and that 'P(> 6) = I - 'P(> 6) thus 
\dP(< 6)/dM\ = \dP(> S)/dM\. Thus the void 
PDF as a function of \6y\ can be obtained from 
the PDF of MVJ (Sll or LMSV [9], provided 
one keeps track of the sign of each term. For 
example in the LMSV approximation the void 
distribution function becomes [|32ll : 
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where the expression is reported as a function 
of the smoothing radius rather than the mass, 
since a void Lagrangian radius is probably eas- 
ier to determine than its mass. 

Note that while a positive skewness (/nl > 0) 
boosts the number of halos at the high mass end 



(and slightly suppress the number of low-mass 
halos), it is a negative skewness that will in- 
crease the voids size distribution at the largest 
voids end (and slightly decrease it for small 
void sizes). Ref. [l32ll concluded that the abun- 
dance of voids is sensitive to non-Gaussianity: 
\6y\ is expected to be smaller than 6c by a factor 
2 to 3. If voids probe the same scales as ha- 
los then they should provide constraints on /nl 
2 to 3 times worse. However voids may probe 
slightly larger scales than halos: in many non- 
Gaussian models, S increases with scales 
(see e.g.. Fig. [B, compensating for the thresh- 
old. 

The approach reviewed here provides a 
rough estimate of the fractional change in abun- 
dance due to primordial non-Gaussianity but 
will not provide reliably the abundance itself. It 
is important to stress here that rigorously quan- 
titative results will need to be calibrated on cos- 
mological simulations and mock survey cata- 
logs. 

4 Effects on the Halo power spectrum 

Recently, Refs. Il43ll46l showed that primordial 
non-Gaussianity affects the clustering of dark 
matter halos (i.e., density extrema) inducing a 
scale-dependent bias for halos on large scales. 
This can be seen for example by considering 
halos as regions where the (smoothed) linear 
density field exceeds a suitable threshold. All 
correlations and peaks considered in the section 
are those of the initial density field (linearly ex- 
trapolated to the present time). Thus for exam- 
ple in the Gaussian case (471 SSI |49|| for high 
peaks we would have the following relation be- 
tween the correlation function of halos of mass 
M, ^hMi^) ^iid that of the dark matter distri- 
bution smoothed on scale R, corresponding to 
mass M, ^^(r): 



(23) 



where bi - 6clcr\ denotes the Lagrangian bias, 
although more refined expressions can be found 
in e.g., L50J and HSU. 



The Lagrangian bias appears here because 
correlations and peaks are those of the initial 
density field (lineary extrapolated). Making the 
standard assumptions that halos move coher- 
ently with the underlying dark matter, one can 
obtain the final Eulerian bias as Be = I + bi, 
using the techniques outlined in ||52l . [|53l . [|50ll 
and [5T]. 

The two-point correlation function of regions 
above a high threshold has been obtained, for 
the general non-Gaussian case, in ll54l . Il55l and 

m-. 

^ft,M(|xi-X2|) = -l+expm (24) 
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where v = ScCTr. For large separations the ex- 
ponential can be expanded to first order. This 
is what we will do in what follows but we will 
comment on this choice below. 

For small non-Gaussianitie^, we can keep 
terms up to the three-point correlation func- 
tion obtaining that the correction to the 
halo correlation function, A^/, due to a non-zero 
three-point function is given by: 

A^/, = ^ [^f (Xi,X2,X2) +^f\xi,Xi,X2)] 



= 44^\xi,Xi,X2) (26) 

For a general bispectrum B(ki ,k2,ks) this yields 
a correction to the power spectrum (see [,46j for 
steps in the derivation): 
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Effectively that is for values of /nl consistent with ob- 
servations 
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Fig. 4 The scale-dependence of the large-scale 
halo bias induced by a non-zero bispectrum for 
different types of non-Gaussianity. The dashed 
line corresponds to the local type and the dot- 
dot-dot-dashed to equilateral type. Figure re- 
produced from [15711 ■ 



where we have made the substitution a. = kj + 
k^ + Ikik/u. Here Mr = WrM. The effect on 
the halo bias is Ab^/b^^ = and thus 



2 P 



&f L = 1 + 



Ac(z) 



[l+6,iz)Mk)] , (28) 



where the expression for /3 can be obtained 
by comparing to Eq. (ITTl) . The term 
A,.(z)/[cr^D^(z)] ^ b^ - I can be recognized 
as the Gaussian Lagrangian halo bias. 

So far the derivation is generic for all types 
of non-Gaussianity specified by a given bispec- 
trum. We can then consider specific cases. In 
particular for local non-Gaussianity we obtain: 
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Thus AbiJbh is 2/nl times a redshift- 
dependent factor Ac(z)/Z)(z) = 6c(z), times a k- 
and mass-dependent factor. The function l3R{k) 
is shown as the dashed line in Fig. dH). This re- 
sult for the local non-Gaussianity has been de- 
rived in at least three other ways. Ref. [i43il 
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generalize the Kaiser [|47ll argument of high 
peak bias for the local non-Gaussianity. Start- 
ing from V^O = + 2/nl [(pV^^ + IV^p] 
where near peaks | V^p is negligible they obtain 
S = (5^NL=o + 2/nl0]. The Poisson equation to 
convert (pin 6 then gives the scale-dependence. 
More details are presented elsewhere in this is- 
sue. 

Ref. f5S\ work in the peak-background 
split. This approach is especially useful to 
understand that it is the coupling between 
very large and small scales introduced by lo- 
cal (squeezed-cofigurations) non-Gaussianity 
to boost (or suppress) the peaks clustering. In 
this approach, the density field can be writ- 
ten as p(x) = p(l + 6i + 6s) where 6i denotes 
long wavelength fluctuations and 6s short wave- 
length fluctuations. 6i is the one responsible 
for modulating halo formation (i.e. to boost 
peaks above the threshold for collapse), so the 
halo number density is n = n{l + bi6i) and 
bL = n~^dn/d6i. 

In the local non-Gaussian case they decom- 
pose the Gaussian field as a combination of 
long and short wavelength fluctuations = 

0, + 0, thus O = 0/ + /NL0f + (1 + 2/nl0/)<^.v + 

fNL((>l + const.. Also in this non-Gaussian case 
one can split the density field 6 in 6i and 6s 
and relate this to /nl (it is easier to work in 
Fourier space): 6i(k) = cc(k)(f>i(k) and 6s = 

a(fc)[(l+2/NL0/)0.+/NL0'] = (X{k)[Xi(f>s+X2(f>^] 

the last equality giving the definition of Xi and 
X2. Note that 6s cannot be ignored here be- 
cause (pi enters in Xi, in other words, local non- 
Gaussianity couples long and short wavelength 
modes. The local halo number density is now 
function of 61, Xi and X2 yielding the following 
result for the Lagrangian halo bias: 
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where (x{k) encloses the scale-dependence of 



the eff'ect. Ref. ||59l rederives the ellipsoidal 
collapse for small deviations from Gaussianity 
of the local type. They find that a non-zero /nl 
modifies the threshold for collapse, the modifi- 
cation is proportional to /nl- This should sound 
familiar from §3. They then use the definition 
hi = n~^dn/d6c keeping track of the fact that 6 
is "modulated" by /nl- 

The efl^ect of the non-Gaussian halo bias on 
the power spectrum is shown in Fig.© where 
the points are measurements from an N-body 
simulations of [33J (see figure caption for more 
details). 

The above result 
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f ~1 

where finik) = ft^L^f^' (k), can be improved in 
several ways. 

First of all, we have not made any distinction 
between the redshift at which the object is being 
observed (z^) and that at which is being formed 
(z/). Except for the rarest events this should be 
accounted for. The Gaussian Lagrangian bias 
expression used so far is an approximation, a 
more accurate expression is |[52l |53l [50ll 
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Then, the halo bias expressions are derived 
within the "classical" Press & Schechter the- 
ory, as we have seen in §3, subsequent im- 
provements on the mass function can be seen 
as a correction to the collapse threshold. In 
the expression for the Gaussian halo bias b1 = 
rT^dff Id6c one can consider mass functions 
that are better fit to simulations than the stan- 
dard Press & Schecter one obtaining: 
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Fig. 5 Effect of the non-Gaussian halo bias on the power spectrum. In the left-top panel we show 
the halo-matter cross-power spectrum for masses above 10^^ Mq at z = 1.02. The left-bottom 
panel shows the ratio of the non-Gaussian to Gaussian bias. Figure reproduced from [33 |. The /nl 
values reported in the figure legend should be interpreted as f^l^. On the right panel we show the 
expected effect and error-bars for the large-scale power spectrum for a survey like LSST. Figure 
reproduced from [|26l . 



The parameters q and p account for non- 
spherical collapse and fit to numerical simula- 
tions yield q ~ 0.75, p = 0.3 e.g., [i34il . In this 
expression the term in the second line is usually 
sub-dominant. The term "-1/5^" in the first 
line is known as "anti-bias", and it becomes 
negligible for old halos Zf «; z^. Note that by 
including the anti-bias correction in bg of Eq. 
|3n one recovers the "recent mergers" approxi- 
mation of Ref. [|58l . 

The same correction should also apply to the 
non-Gaussian correction to the halo bias: 

Ab = f^i^q'6,(b'^ - l)/3f'=\k) . (34) 

where q' should coincide with q above; it 
can be calibrated to N-body simulations and is 
found indeed to be q' = 0.75 [|33l . The non- 
Gaussian halo bias prediction and results from 
N-body simulations with local non-Gaussianity 
are shown in Fig. Q. 

Finally one may note that for /nl large and 
negative, Eq. (|27l[28l) would formally yield bf'^ 



and Ph(k) negative on large enough scales. This 
is a manifestation of the breakdown of the ap- 
proximations made: a) all correlations of higher 
order than the bispectrum were neglected: for 
large NG this truncation may not hold; b) The 
exponential in Eq. (I24l) was expanded to linear 
order. This however could be easily corrected 
for, remembering that the P(k) obtained from 
Eq. (ITT] ) is in reality the Fourier transform of 
X, the argument of the exponential. One would 
then compute the halo correlation function us- 
ing Eq. (I24l) and Fourier transforming back to 
obtain the halo power- spectrum. 

So far we have concentrated on local non- 
Gaussianity, but the expression of Eqs. (I27l[28]) 
is more general. Using this formulation, Ref. 
[|57l computed the quantity (i^^^'^ik) for sev- 
eral types of non-Gaussianity (equilateral, lo- 
cal and enfolded); this is shown in Fig. (U). 
It is clear that the non-Gaussian halo-bias ef- 
fect has some sensitivity to the bispectrum 
shape, for example the effect for the equilateral 




Fig. 6 The quantity ^b/b as function of k, for simulation snapshots at z = 0.44, 1.02 and 1.54. 
Simulation outputs and theory lines are shown for f^i = +100 and f^t = ±200. Figure reproduced 
from Il33l . The /nl values reported in this figure legend should be interpreted as f^l^. 



type of non-Gaussianity is suppressed by orders 
of magnitude compared to the local-type and 
the flattened case is somewhere in the middle. 
Fig. @) also shows a type of non-Gaussianity 
arising from General-relativistic (RG) correc- 
tions on scales comparable to the Hubble ra- 
dius. Note that perturbations on super- Hubble 
scales are initially needed in order to "feed" 
the GR correction terms. In this respect the 
significance of this contribution is analogous 
to the well-known large-scale anti-correlation 
between CMB temperature and E-mode po- 
larization: it is a consequence of the proper- 
ties of the inflationary mechanism to lay down 
the primordial perturbations. This effect has 
the same magnitude as a local non-Gaussianity 
with/NL > 1. 

The next logical step is then to ask how well 
present or forthcoming data could constrain 
non-Gaussianity using the halo-bias effect. It is 
interesting to note that surveys that aim at mea- 
suring Baryon Acoustic Oscillations (BAO) in 
the galaxy distribution to constrain dark energy 
are well suited to also probe non-Gaussianity: 
they cover large volumes and their galaxy num- 
ber density is well suited so that on the scale 
of interest (both for BAO and non-Gaussianity) 
shot noise does not dominate the signal. Pho- 



tometric surveys are also well suited: as the 
non-Gaussian signal is localized at very large 
scales and is a smooth function of k, the photo- 
z smearing effects are unimportant. 

The theory developed so far describes the 
clustering of halos while we observe galaxies. 
Different galaxy populations occupy dark mat- 
ter halos following different prescriptions. If 
we think in the halo-model framework (e.g., 
[f60l and references therein) at very large scales 
only the "two halo" contribution matters and 
the details of the halo occupation distribution 
(the "one-halo" term) become unimportant. 

What is important to keep in mind is that 
the effect of the non-Gaussianity parameter one 
wants to measure, /nl, is fully degenerate with 
the value of the Gaussian (small scales) halo 
bias. Fig. ^ shows the dependence of the non- 
Gaussian correction on the Gaussian bias. 

Thus highly biased tracers will show a a 
larger non-Gaussian effect for the same /nl 
value. Of course for a given cosmological 
model the Gaussian bias can be measured accu- 
rately by comparing the predicted dark matter 
power spectrum with the observed one. Alter- 
natively, two differently biased tracers can be 
used in tandem to disentangle the two effects 
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Fig. 7 Non-Gaussian halo bias correction as 
function of the gaussian halo bias. Figure re- 
produced from |[33l .The /ml values reported in 
the figure legend should be interpreted as f^^^ . 

Since clustering amplitude may depend on 
the entire halo history, it becomes then inter- 
esting to model in details the dependence of the 
effect on the halo merger tree (Reid et al., in 
preparation). 

4.1 Outlook for the future How well can 
this method do to constrain primordial non- 
Gaussianity compared with the other tech- 
niques presented here? The Integrated Sachs 
Wolfe (ISW) effect offers a window to probe 
clustering on the largest scales (where the sig- 
nal is large); on the other hand, a measure- 
ment of clustering of tracers of dark matter ha- 
los is a very direct window into this effect. A 
Fisher matrix approach [63l|59l[58l shows that 
the ISW signal is weighted at relatively low 
redshift (where dark energy starts dominating) 
while the non-Gaussian signal grows with red- 
shift, thus making the shape of the halo power 
spectrum a more promising tool. An overview 
of current constraints from different approaches 
can be found in table [T] and future forecasts in 
table [2l for non-Gaussianity of the local type. 
Large, mass-selected cluster samples as pro- 
duced by SZ-based experiments will provide a 
optimally suited data-set for this technique (see 



e.g. HMD. 

While for a given /nl model such as the lo- 
cal one, methods that exploit the non-Gaussian 
bias seem to yield the smallest error-bars for 
large-scale structure, it should be kept in mind 
that the bispectrum can be used to investigate 
the full configuration dependence of /nl and 
thus is a very powerful tool to discriminate be- 
tween different type of non-Gaussianity. In 
addition CMB -bispectrum and halo bias test 
non-Gaussianity on very large scales while the 
large scale structure bispectrum mostly probes 
mildly non-linear scales. As primordial non- 
Gaussianity may be scale-dependent, all these 
techniques are highly complementary. 

The above estimates assume that the under- 
lying cosmological model is known. The large- 
scale shape of the power spectrum can be af- 
fected by cosmology. Carbone et al. (in prep.) 
explore possible degeneracies between /nl and 
cosmological parameters. They find that the pa- 
rameters that are most strongly correlated with 
/nl are parameters describing dark energy clus- 
tering, neutrino mass and running of the pri- 
mordial power spectrum spectral slope. For 
surveys that cover a broad redshift range the 
error on /nl degrade little when marginaliz- 
ing over these extra parameters: the peculiar 
redshift dependence of the non-Gaussian signal 
lifts the degeneracy. 

5 Conclusions 

A natural question to ask at this point may be 
"what observable will have better chances to 
constrain primordial non-Gaussianity?" 

In principle the abundance of rare events is a 
very powerful probe of non-Gaussianity; how- 
ever, in practice, it is limited by the practical 
difficulty of determining the mass of the ob- 
served objects and its corresponding large un- 
certainty in the determination. This point is 
stressed e.g., in O. With the advent of high- 
precision measurements of gravitational lens- 
ing by massive clusters, the mass uncertainty, at 
least for small to moderate size clusters samples 
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Table 1 Current recent 2 
Data/method 



cr constraints on local fnl 

/nl 



reference 



Photometric LRG - bias 
Spectroscopic LRG- bias 

QSO - bias 

combined 



z:o+54+101 

7A+74+I39 
' -83-191 
o +26+47 
°-37-77 
90+23+42 
-24-57 



Slosar et al. 2008 
Slosar et al. 2008 
Slosar et al. 2008 
Slosar et al. 2008 



1 +647+755 
^^-^-337-1157 

236+ 127(2 -(t) 



NVSS-ISW 
NVSS-ISW 



Slosar et al. 2008 
Afshordi&ToUey 2008 



30 + 84 Spergel et al. (WMAP) 2007 

32 + 68 Creminelli et al. 2007 

87 + 60 Yadav & Wandelt 2008 

38 + 42 Smith et al. 2009 

5 1 + 60 Komatsu et al. (WMAP) 2008 

-57 + 121 Komatsu et al. (WMAP) 2008 



WMAP3-Bispectrum 
WMAP3-Bispectrum 
WMAP3-Bispectrum 
WMAP-Bispectrum 
WMAP5-Bispectrum 
WMAP5-Minkowski 



Table 2 Forecasts 1 - cr constraints on local /nl 



Data/method 


A/NL(l-cr) 


reference 


BOSS-bias 


18 


Carbone et al. 2008 


ADEPT/Euclid-bias 


1.5 


Carbone et al. 2008 


PANNStarrs -bias 


3.5 


Carbone et al. 2008 


LSST-bias 


0.7 


Carbone et al. 2008 


LSST-ISW 


7 


Afshordi& ToUey 2008 


BOSS-bispectrum 


35 


Sefusatti & Komatsu 2008 


ADEPT/Euclid -bi spectrum 


3.6 


Sefusatti & Komatsu 2008 


Planck-Bispectrum 


3 


Yadav et al . 2007 


BPOL-Bispectrum 


2 


Yadav et al . 2007 



can be greatly reduced. Forthcoming Sunyaev- 
Zeldovich experiments will provide large sam- 
ples of mass-selected clusters which could then 
be followed up by lensing mass measurements 
(see e.g., [[65l [66l|). So far there is only one 
very high redshift (z = 1 .4) very massive M ^ 
8 X 10^"^ Mq with high-precision mass determi- 
nation via gravitational lensing {WT\. Ref. [168II 
pointed out that this object is extremely rare, 
for Gaussian initial conditions there should be 
0.002 such objects or less in the surveyed area, 
which is uncomfortably low probability. But 
the cluster mass is very well determined: a non- 
Gaussianity still compatible with CMB con- 
straints could bring the probability of observing 
of the object to more comfortable values. This 



result should be interpreted as a "proof of prin- 
ciple" showing that this a potentially powerful 
avenue to pursue. 

The measurement of the three-point corre- 
lation function allows one to map directly the 
shape-dependence of the bispectrum. For large- 
scale structures the limiting factors are the large 
non-Gaussian contribution induced by gravita- 
tional evolution and the uncertainly of the non- 
linear behavior of galaxy bias. 

The halo-bias approach can yield highly 
competitive constraints, but it is less sensitive to 
the bispectrum shape. Still, the big difference in 
the magnitude and shape of the scale-dependent 
biasing factor between different non-Gaussian 
models implies that the halo bias can become a 
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Table 3 Forecasted non-Gaussianity con- 
straints: A) [20] B) [63] C) den [Toll E) (53) 
e.g., [15] 



CMB Bispectrum Halo bias 


type NG 


Planck 


BPol Euclid 


LSST 


l-cr errors 


Local 


3A) 


2A) 155) 


0.7^) 


Equilateral 


25^) 


14C) 




Enfolded 


(910 


(910 39^) 


18^) 


#cr detection 


GR 


N/A 


N/A 1^) 


2£) 


secondaries 


3F) 


5^) N/A 


N/A 



useful tool to study shapes when combined with 
e.g. measurements of the CMB bispectrum. Ta- 
ble [3] highlights this complementarity. For ex- 
ample, one could envision different scenarios. 

If non-Gaussianity is local with negative f^i 
and CMB obtains a detection, then the halo bias 
approach should also give a high- significance 
detection (GR correction and primordial con- 
tributions add up), while if it is local but with 
positive f^L, the halo-bias approach could give 
a lower statistical significance for small f^i as 
the GR correction contribution has the opposite 
sign. 

If CMB detects /nl at the level of ~ 10 and of 
a form that is close to local, but halo bias does 
not detect it, then the CMB bispectrum is given 
by secondary effects. 

If CMB detects non-Gaussianity but is not of 
the local type, then halo bias can help discrim- 
inate between equilateral and enfolded shapes: 
if halo bias sees a signal it indicates enfolded 
type,; if halo bias does not see a signal it in- 
dicates equilateral type. Thus even a non- 
detection of the halo-bias effect, in combination 
with CMB constraints can have an important 
discriminatory power. 

In any case, if the simplest inflationary sce- 
nario holds, for surveys like Euclid and LSST, 
the halo-bias approach is expected to detect a 
non-Gaussian signal very similar to the local 



type signal with an amplitude of f^i ~ -1.5 
which is due to large- scales GR corrections to 
the Poisson equation. This effect should leave 
no imprint in the CMB: once again the com- 
bination of the two observable can help enor- 
mously to discriminate among models for the 
origin of cosmological structures. 

In addition we should bear in mind that non- 
Gaussianity may be scale-dependent. In fact for 
models like DBI inflation it is expected to be 
scale-dependent. A proposed parameterization 
of the scale-dependence of non-Gaussianity is 
given by : 

— j F(^i,4,4) (35) 

where K = (kikjhy'^ JTOl, kp denotes the 
pivot and n„g the slope or running of non- 
Gaussianity , although other authors prefer to 
use K = (ki+k2 + h)/?) [|9l|59l as for squeezed 
configurations K 0. It is still an open is- 
sue which parameterization is better in practice. 

In any case different observables probe dif- 
ferent scales (see Fig 18]) and their complemen- 
tary means that "the combination is more than 
the sum of the parts". 

What is clear, however, is that the thorny sys- 
tematic effects that enter in all these approaches 
will require that a variety of complementary 
avenues be taken to establish a robust detection 
of primordial non-Gaussianity. 
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Fig. 8 Scale-dependent /nl and scales probed by different approaches mentioned here. The solid 
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